Power spectrum for fractal distributions 
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We study the behaviour of the power spectrum (PS) in the case of fractal structures. We 
show that in this case the main observational features of the PS, the large scale flattening and 
the scaling of the amplitude with sample depth, are related to finite size effects, due to the fractal 
nature of galaxy distribution in the sample. Comparing with the recent results of the PS for the 
CfA2 we conclude that this catalog is consistent with the fractal description. 

Identifying the scale at which our Universe becomes homogeneous, if any, is a 
crucial task of contemporary cosmology, and a very debated one, especially for what 
concerns the homogeneity of the luminous matter at the present. Essentially all the 
currently elaborated models of galaxy formation assume large scale homogeneity 
and predict that the galaxy power spectrum (PS) decreases both toward small 
scales and toward large scales, with a turnaround at some scale A/ that can be 
taken as separating "small" from "large" scales. Then, because of the assumption 
of homogeneity, the power spectrum amplitude should be independent of the survey 
scale, any residual variation being attributed to luminosity bias (or to the fact that 
the survey scale has not yet reached the homogeneity scale). However, the crucial 
clue to this picture, the firm determination of the scale A/, is still missing, although 
some surveys do indeed produce a turnaround scale around 100 h~ 1 Mpc (Baugh & 
Efstathiou 1992; Feldman et aL1994). Recently, the CfA2 survey analyzed by Park 
et qZ.(1994; PVGH) (and confirmed by SSRS2 - Da Costa et o/.(1994, DVGHP)), 
showed a n = —2 slope up to ~ ZQh~ 1 Mpc 1 a milder n m — 1 slope up to 200 
/i _1 Mpc, and some tentative indication of flattening on even larger scales. They also 
find that deeper subsamples have higher power amplitude, i.e. that the amplitude 
scales with the sample depth. 

In this paper we argue that both features, bending and scaling, are a manifes- 
tation of the finiteness of the survey volume, and that they cannot be interpreted 
as the convergence to homogeneity, nor to a power spectrum flattening. The sys- 



tematic effect of the survey finite size is in fact to suppress power at large scale, 
mimicking a real flattening. We show that even a fractal distribution of matter, 
i.e. a distribution which never reaches homogeneity, shows a sharp flattening in the 
PS (even when the correction proposed by Peacock & Nicholson (1991) is applied 
to the data), and that its amplitude depends on the survey size (see Sylos Labini 
& Amendola, 1995). The reason of such behavior is that the standard power spec- 
trum (hereafter SPS) measures directly the contributions of different scales to the 
galaxy density contrast Sp/p. It is clear then that the density contrast, and all the 
quantities based on it, is meaningful only when one can define a constant density, 
i.e. reliably identify the sample density with the average density of all the Universe. 
When this is not true, and we argue that is indeed an incorrect assumption in all 
the cases investigated so far, a false interpretation of the results may occur, since 
both the shape and the amplitude of the power spectrum depend on the survey size. 

Let us recall the basic notation of the power spectrum analysis. Following 
Peebles (1980) we imagine that the Universe is periodic in a volume V u , with V u 
much larger than the (assumed) maximum correlation length. The survey volume 
V 6 V u contains N galaxies at positions fl, and the galaxy density contrast is 
S(r) = [n(r)/h] — 1 where it is assumed that exists a well defined constant density 
n, obtained averaging over a sufficiently large scale. The density function can be 
described by a sum of delta functions: n(r) = J2iLi S^ 3 \r ~ • Expanding the 
density contrast in its Fourier components we have 

5 % = 1/Nj2^ rl ~W(k), 1 

jeV 

where W(k) = V^ 1 J drW (r)e tkr is the Fourier transform of the survey window 
W(r), defined to be unity inside the survey region, and zero outside. If £(r) is 
the correlation function of the galaxies, (£(r) =< n(r)n(0) > /n 2 — 1) the true PS 
P{k) is defined as the Fourier conjugate of the correlation function £(r). Because 
of isotropy the PS can be simplified to 

P(k) =4tt J t(r)sin(kr)/(kr)r 2 dr . 2 

The variance of <% is (Peebles 1980) < |<%| 2 >= A^ 1 + V _1 P(fc) . i.e. the sum of a 
shot noise term and of the true PS convolved with a window function (Sylos Labini 
& Amendola, 1995) 

We apply now this standard analysis to a fractal distribution. Consider a self- 
similar system, where the number of points inside a certain radius r scales according 
to the mass-length relation (Mandelbrot, 1982) N(r) = Br D , with D < 3 (the case 
D = 3 corresponds to the homogenous distribution). Then we have the following 
correlation function 

^(r) = [(3- 7 )/3](r/i2 a )-T-l, 3 

where 7 = 3 — D, On scales larger that R s the £(r) cannot be calculated without 
making assumptions on the distribution outside the sampling volume (in particular, 
assumption of homogeneity), which is just what we want to avoid. When the survey 
volume is not spherical, the scale R s is of the order of the largest sphere completely 
contained inside the survey. Both the amplitude and the shape of £(r) are therefore 



scale-dependent in the case of a fractal distribution (CP92). inside a sphere of 
radius R s turns out to be 

P(k) = (a(k, R s )R 3 ~ D )/(k D ) - (b(k, R s ))/(k 3 ) . 4 

Notice that the integral has to be evaluated inside R s because we want to com- 
pare P(k) with its estimation in a finite size spherical survey of scale R s . In the 
general case, we must deconvolve the window contribution from P(k); R s is then a 
characteristic window scale. The previous equation shows the two scale-dependent 
features of the PS. First, the amplitude of the PS depends on the sample depth. 
Secondly, the shape of the PS is characterized by two scaling regimes: the first one, 
at high wavenumbers, is related to the real fractal dimension, while the second one 
arises only because of the finitcness of the sample. In the case of D = 2 one has: 
a = (47r)/(3)(2 + cos(fci? s )) , and b = A-k sa\{kR s ) . The PS is then a power-law 
with exponent —2 at high wavenumbers, it flattens at low wavenumbers and reaches 
a maximum at k w 4.3/ R s , i.e. at a scale A « 1.45i? s . In a real survey, things are 
complicated by the window function, so that the flattening (and the turnaround) 
scale can only be determined numerically. 

To avoid the mean density normalization, which gives misleading results in 
fractals distributions, we introduce now the scale-independent PS (SIPS) of the 
density p(r), a quantity which gives an unambiguous information of the statisti- 
cal properties of the system. We first introduce the density correlation function 
G(r) =< p(x + r)p(x) >= Ar~^~ D ^ , where the last equality holds in the case of 
a fractal distribution with dimension D, and where A is a constant Defining the 
SIPS as the Fourier conjugate of the correlation function G(r), one obtains that in 
a finite spherical volume II(fc) ~ A'k~ D , (where A' = — cos(kR s )) if D = 2) 
so that the SIPS is a single power law extending all over the system size, without 
amplitude scaling (except for kR s < 1). In analogy to the procedure above, we 
consider the Fourier transform of the density p^ = V^ 1 Yljev e~ tkXj , and its vari- 
ance < \p^\ 2 >= V _1 II(fc) + N^ 1 , where H(k) is the same of the previous case but 
with < \pg\ 2 > instead of < \5g | 2 >. 

In Fig. 1 we show the observational results of PVGH and compare them with the 
PS of a fractal distribution with D = 2, adopting the same technique used for real 
redshift surveys (Peacock & Nicholson 1991; Fisher et al. 1993; PVGH; DVGHP). 
Let us summarize the observational results of PVGH, by confronting them with 
the PS for a fractal distribution (Fig. 1): i) for k > 0.25 (A < 25/1' 1 Mpc) the 
PS in a volume limited sample is very close to a power law with slope n = —2.1. 
In our view, this is the behaviour at high wavenumbers connected with the real 
fractal dimension, ii) For 0.05 < k < 0.2 {120h~ 1 Mpc > A > SOh^Mpc) and 
the spectrum is less steep, with a slope about —1.1. This bending is, in our view, 
solely due to the finite size of the sample, iii) The amplitude of the volume limited 
subsample CfA2-130 PS is ~ 40% larger than for CfA2-101. This linear scaling of 
the amplitude can be understood again considering that the sample is fractal with 
D = 2. The same behaviour has been found in the analysis of the £( r ) : ro scales 
linearly with the sample depth according to Eq.3 (CP92). 

The authors (PVGH) explain this fact considering the dependence of galaxy 
clustering on luminosity: brighter galaxies correlate more than fainter ones. It 
is certainly possible that both mechanisms, the luminosity segregation and the 



FIGURE 1. Comparison of power spectra of fractal distribution (triangles) with 
the CfA2 survey (squares). In the top panel, we plot the PS of the subsamplc 
CfA2-130 (PVGH) along with the PS of our artificial fractal distribution (without 
the error bars for clarity). In the bottom panel, we plot CfA2-101 (PVGH) and a 
subsamplc of the same fractal as above, with a correspondingly scaled depth. 

intrinsic self-similarity of the distribution, are correct, and each one explains part 
of the scaling. However, PVGH do not detect such a luminosity segregation for 
the two largest subsamples, CfAlOl and CfA130, to which we are comparing our 
analysis here. It seems therefore that the amplitude scaling at these scales can be 
entirely attributed to the fractal scaling. Our conclusion is then that the fractal 
nature of the galaxy distribution can explain, to the scales surveyed so far, the shift 
of the amplitude with sample depth of the PS and of £ (r) . 

Finally we stress that the fractal dimension of the galaxy clustering rises from 
D <~ 1.4 for CfAl to D ~ 2 for CfA2, in agreement with the result of other indepen- 
dent surveys: Perseus-Pisces (Guzzo et a/.1992; Sylos Labini et al.1995), and ESP 
(Pietronero & Sylos Labini 1994, Baryshev et al., 1994). 
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